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The decays — >■ Ds*'"'"Ds*'~ are reconstructed in a data sample corresponding to an integrated 
himinosity of 6.8 fb~^ collected by the CDF II detector at the Tevatron pp collider. All decay 
modes are observed with a significance of more than 10 a, and we measure the production rate 
times Bg _Ds*'^-Ds*' branching ratios relative to the normalization mode B'' Df D~ to be 
0.183±0.021±0.017 for B° D+D', 0.424 ±0.046 ±0.035 for B° D'^^Df, 0.654±0.072±0.065 
for B° D*+D*-, and 1.261 ± 0.095 ± 0.112 for the inclusive decay B'^ where the 

uncertainties are statistical and systematic. These results are the most precise single measurements 
to date and provide important constraints for indirect searches for non-standard model physics in 
Bs mixing. 
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A meson can oscillate into its antiparticle via 
second order weak interaction transitions, which make 
its time evolution sensitive to contributions from new 
physics processes. Such contributions are not well con- 
strained yet and might be responsible for the deviation 
from the standard model reported in Ref. The 
cigcnstates with defined mass and lifetime, and -6°^, 
are linear combinations of the 1?° and states and, in 
the standard model, correspond in good approximation 
to the even and odd CP eigenstates, respectively. In 
the absence of substantial CP violation, a sizable decay 
width difference between the light and heavy mass eigen- 
states, AFs = TsL — TsH, arises from the fact that decays 
to final states of definite CP are only accessible by one 
of the mass eigenstates. The dominant contribution to 
AFs is believed to come from the ^ d':*^'^ d':*^ de- 



Lubbock, TX 79609, USA, '^'^Universidad Tecnica Federico Santa 
Maria, llOv Valparaiso, Chile, '^''Yarmouk University, Irbid 211- 
63, Jordan, '""'University of Warwick, Coventry CV4 7AL, United 
Kingdom, 



4 



cays 0, which are predominantly CP-even and saturate 
ATs under certain theoretical assumptions 0] , result- 
ing in the relation 



Ar, 



r, + Ar,/2' 



(1) 



where Ts = {TsL + Tsh)/'^ 01 • However, three-body 
modes may provide a significant contribution to AF^ [6| . 

A finite value of AF^ improves the experimental sen- 
sitivity to CP violation because it allows one to dis- 
tinguish the two mass eigenstates via their decay time 
distribution. Furthermore, the — > Z?s*''^_Di*-' de- 
cays could be used in future to measure directly the 
lifetime of the CP-even eigenstate, which would com- 
plement the CP-odd eigenstate lifetime measurement in 
^ J/V'/o(980) decays and provide additional in- 
formation in the search for new physics contributions to 
CP violation in the -B° system. 

The — > Ds*''^£'i*^ decay modes have been previ- 
ously studied by the ALEPH, CDF, DO, and Belle col- 
laborations d - 11 1 . The current world average branching 
ratios [l^ , which do not yet include the latest preliminary 
Belle results [H, are 6(5° ^ D+Dj) = (1.04j:[;-^9) %, 
B{B°^ D^'^Df) = (2.8 ± 1.0) %, B{B° D*+D*-) = 
(3.1 ± 1.4) %, and B{B° 7?i*'+i?i*'^) = (4.5 ± 1.4) %. 

In a data sample corresponding to an integrated lumi- 
nosity of 6.8 fb~^ recorded by the CDF II detector at the 
Tevatron pp collider we reconstruct B'^ — !• ui*"^^ ui*^ 
decays with A'+A'^7r+. For the first time in 

this channel, the acceptance is calculated using a Df 
Dalitz model instead of a simple two-body decay model. 
The photon and the neutral pion from the D*^ — !■ Dfj 
and Z?*"*" — >■ Dfir'^ decays are not reconstructed because 
of their low detection efficiency. In a simultaneous fit 
to the reconstructed S^^^ meson invariant mass spec- 
tra we measure the 5° production rate times 5° — ?> 
Di*^+Di*^- bi 
tion mode B° D+D- 



Ds Ds branching ratios relative to the normaliza- 



X 



fs B{B^, ^ X) 
fdB{B° ^ D+D-)' 



(2) 



for X = DtDj, D*+Df, D*+ D*' , and the inclusive 
where fs/fd is the relative rate of produced 
mesons. 



B° to B° 

The components of the CDF II detector 1^ most rel- 
evant for this analysis are the tracking systems located 
inside a solenoid that provides a 1.4 T magnetic field. 
Charged particles' trajectories (tracks) are reconstructed 
in layers of silicon-strip sensors located between radii of 
1.5 cm and 28 cm from the beam line and an open-cell 
drift chamber (COT) with a radial extension from 40 to 
137 cm. Tracks with a pseudorapidity \ri\ < 1.0 pass 
the full radial extent of the COT. Kaons and pions are 
statistically identified by measurements of the ionization 



energy loss in the COT and information from the time-of- 
fiight system located between the COT and the solenoid. 
The events for this analysis are selected online by identi- 
fying pairs of tracks detected in the COT and the silicon 
detector |15| . Minimal requirements on the momenta and 
the displacement of the tracks and the reconstructed de- 
cay vertex from the primary vertex are imposed. 

We reconstruct D+ — >■ K+K~'k+ and D+ — >■ i4r~7r+7r+ 
decays from combinations of three tracks with appropri- 
ate charge and mass hypothesis assignments, fitted to a 
common vertex. Because the Df — !• K+K~Tr+ decay 
proceeds mainly via 07r+ and K*'^K+, we select can- 
didates with 1.005 < m{K+K-) < 1.035 GeV/c^ and 
0.837 < m{K-Tr+) < 0.947 GcV/c^, centered on the 
known (j> and K*'^ masses, respectively. According to the 
Df — 7> K+K~TT+ Dalitz structure [16j this requirement 
has a signal acceptance of about 75 % while covering only 
14 % of the phase space and thus increasing the signal- 
to-background ratio. In the following we will denote the 
selected K+K~ and K~7t+ combinations as cp and K*'^, 
respectively, since the dominant contributions come from 
these resonances. However, we implicitly include con- 
tributions from other resonances and interference effects 
when using these terms. 

Pairs of D+ — ^ (j)7T+ or Df K*^K+ candidates and 
Dj (j>Tr~ candidates are combined to form B^ can- 
didates and fitted to a common vertex. Combinations 
where both charm mesons decay into a K*'^ mode are not 
considered because of the low signal-to-background ra- 
tio. Candidate B° mesons are reconstructed from D+D~ 
combinations where both Df decay modes are used. 

To reject background-like events, requirements are 
placed on track quality variables, B meson momentum, 
reconstructed D meson masses, vertex fit qualities, and 
vertex displacement significances. To further increase the 
signal purity, two artificial neural networks are used, one 
for candidates with a K*'^ and one for candidates with- 
out. To minimize the systematic uncertainty of the rela- 
tive selection efficiency, the same networks are applied to 
and B° candidates, and only information from the Df 
that is common to both B meson decays is used. The net- 
works are trained on simulated signal events, described 
below, and on background events from the 5.45 to 6.5 
GeV/c^ B mass sideband. The input variables contain 
kinematic, lifetime, fit quality, and particle identification 
information. The B vertex displacement significance in 
the transverse plane gives the largest contribution to the 
discrimination power of both networks. The selection cri- 
teria on the network outputs are chosen such that they 
maximize the significance emc / V-^data, where cmc is the 
-Bj selection efficiency determined from simulation and 
-/Vdata is the number of data events in the B^ signal win- 
dow from 5.343 to 5.397 GeV/c^. 

About 6% of the selected D+{-^ (pTT+)D- can- 

didates also fulfill the B^ selection requirements, where 
the assignment of a D~ daughter track is swapped from 
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pion to kaon. To avoid having the same event entering 
the fit multiple times, we reject each event that is re- 
constructed as -B° candidates from the sample. The 
cross-populations between the two 5° modes and be- 
tween the two i?° modes, respectively, are negligible. The 
selected sample contains about 750 signal events. 

Simulated events are used to determine the reconstruc- 
tion and selection cfhciency. The mesons arc gen- 
erated according to the momentum spectrum measured 
in exclusive B decays and decayed to the considered fi- 

For the B° 
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nal states with the evtgen package 
meson wc assign the lifetime of the eigenstate [Ij 
that coincides with the CF-even eigenstate in the stan- 
dard model. For all the other long-lived charm and bot- 
tom mesons, the world average mean lifetimes 12] are 
used. The i?° — > D*+D*~ decay is a transition of a 
pseudoscalar to two vector mesons and its angular dis- 
tribution is described by three polarization amplitudes. 
Since these amplitudes are unknown, wc take the same 
longitudinal polarization as measured in — D*^D*~ 
decays and a vanishing CP-odd component as default 
values. The world average value [l^l is used for the ratio 
of D*+ D+j to D*+ D+Tr° decays. The dynamics 
of the decay Df — > K^K^tt^ is simulated according to 
the Dalitz structure measured by CLEO . The gen- 
erated events are processed by a GEANt3 based detector 
simulation [l^ and the same reconstruction program as 
applied to real data events. 

The relative branching ratios times produc- 
tion rate are determined in a simultaneous ex- 
tended unbinned maximum-likelihood fit to 



the 
and 

By 

nor- 



{K*'^K~^){K^TT^TT~) invariant mass distributions, 
simultaneously fitting all four distributions, the 
malization of the i?° reflections in the {K*°K'^){(f>Tr^) 
spectrum is constrained by the yields in the high- 
statistics (07r+)(A'+7r~7r^) sample. The components 
of the fit function for each invariant mass distribution 
are fully and partially reconstructed signals, reflections, 
and background. The fully reconstructed i?° and -B" 
signals are parametrized by the sum of two Gaussians 
with relative normalizations and widths derived from 
simulation. To account for discrepancies between data 
and simulation, a factor is introduced for the B^ and 
B^ signal shapes, respectively, that scales the widths 
of the Gaussians and that is allowed to float in the 
fit. The shapes of partially reconstructed signal events 



-){K+TT- 



){K*^K ) arc determined 



and of reflections from 
reconstructed as B^ — > 
from simulation using empirical models. Background 
from random combinations of tracks and other B decays 
is described by an exponential plus a constant function 
with all parameters floated in the flt. 

The yield of fully reconstructed B^ mesons in the fi- 
nal state i, (07r+)(if+7r-7r-) or {K*°K+){K+'k-'!t-), is 



given by 



N^,B{B' ^ DtD-)B{Dt ^ K+K-^+) 



■B{D+ K-7r+Tr+)ef 



(3) 



where N^l is the total number of produced B^ mesons 
and is a free parameter in the fit, the branching ratios 
are taken from Ref. and the efficiency ef is deter- 
mined from simulation. Equivalent expressions are used 
for the yields of partially reconstructed B^ decays with 
an additional branching ratio factor for the D*~^ and D*'^ 
decays. The normalizations of refiections are calculated 
in the same way, but with the efficiencies replaced by 
the mis-reconstruction fractions determined from simu- 
lation. The number of fully reconstructed i?^ mesons in 
the final state i, (07r+)(07r~) or {K*'^K^){<j>Tr~), where 
the decays in the same mode as the from the B° 
decay is given by 



BjPt K+K-TT+) e' 
B{D+ K-TT+TT+) e 



(4) 



with fo^Ds as a free parameter and ^ given by Eq. 
(|3|). Equivalent equations hold for partially reconstructed 
B^ decays. 

Projections of the fit result arc compared to the dis- 
tribution of data events in Fig. [TJ The statistical sig- 
nificance of each signal exceeds 10 a as estimated from 
a likelihood ratio of the fit with and without the signal 
component. 

Systematic uncertainties on the fitted signal yields 
arise from the signal and background models. Because 
the width scale factors of the fully reconstructed signal 
components are allowed to float in the fit, the systematic 
uncertainties of these components are already included 
in the statistical errors. To estimate the systematic ef- 
fect due to the fixed shapes of the partially reconstructed 
signal components and reflections, we repeat the fit mul- 
tiple times with shape parameters randomly varied ac- 
cording to the covariance matrix of the fits of the shapes 
to simulated data. The mean deviations with respect to 
the central values are assigned as systematic uncertain- 
ties. The systematic uncertainties due to the background 
mass model are estimated from the changes in the results 
caused by using a second order polynomial instead of the 
sum of an exponential and a constant function. By ap- 
plying the selection optimization procedure on the nor- 
malization instead of the signal mode we verified that a 
possible selection bias is negligible. 

Systematic effects in the relative efficiency determina- 
tion can be caused by a simulation that does not describe 
the data accurately. One source of systematic uncertain- 
ties is the trigger simulation, which can lead to a discrep- 
ancy in the B meson momentum spectrum. Although 
this effect cancels to first order in the ratio measure- 
ment, it is accounted for by a reweighting of the simulated 
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FIG. 1: Invariant mass distribution of (a) 
D+{ct>-K+)D~{<t)'K-), (b) Bl ^ D+{K-°K+)D7[<t)'K-), 
(c) B" Dt{(l)TT+)D-{K+Tr-Tv-), and (d) B° -> 

TT n ) candidates with the simultane- 
ous fit projection overlaid. The broader structures stem from 
decays where the photon or tt*^ from the O^^gj decay is not 
reconstructed. Misreconstructed signal events in (c) show up 
as reflections in (b). 



events. The systematic uncertainties due to the detector 
simulation are estimated by the shift of the results with 
respect to the case in which this reweighting is not ap- 
plied. The uncertainties on the world average , 
and Df lifetimes are propagated by varying the lifetimes 
in the simulation. For the B2 lifetime, we consider two 



Source 
Signal model 
Background model 
Detector simulation 
B, D lifetimes 
Dalitz model 
Helicity model 
Branching fractions 



fo^D^, /djDs /djdj fr 



(*) 



0.003 
0.001 
0.001 

+0.001 
-0.002 

0.011 
0.001 
0.013 



0.007 
0.004 
0.003 

+0.002 
-0.004 

0.024 
0.005 
0.024 



0.009 
0.030 
0.010 

+0.003 
-0.006 

0.038 
0.012 
0.039 



0.019 
0.033 
0.005 

+0.006 
-0.012 

0.073 
0.008 
0.074 



Total 



0.017 0.035 0.065 0.112 



TABLE I: Overview of systematic uncertainties on the mea- 
sured ratios of branching fractions. 



cases, the la lower bound of the world average short-lived 
eigenstate lifetime and the la upper bound of the mean 

lifetime. The effects on the acceptance induced by 
variations of the — >■ K~^K~Tr~^ Dalitz structure are 
considered by generating different Dalitz model scenar- 
ios, with Dalitz model parameter values varied according 
to the systematic and correlated statistical uncertainties 
of the CLEO Dalitz fit. The uncertainties of the 
Dalitz model have a negligible effect on the result. For 

— > D*^Dl^ decays we investigate the effects of both 
a longitudinal polarization fraction deviating from our 
nominal assumption and a non-zero fraction of the CP- 
odd component fcp- ■ The fraction is varied in the 
simulation according to the uncertainty of the /l mea- 
surement in B*^ D*^D*^ decays [18[. A variation of 
fcp- shows no effect on the 5° — > D*~^D*~ mass line 
shape, fit quality, or measured branching fraction ratios. 
The effect of self cross-feed due to a wrong assignment of 
kaon and pion masses is negligible. 

Further systematic uncertainties arise from external in- 
put quantities. The uncertainties of intermediate and 
final state branching fractions, B{Df /v+/v~7r+), 
B{D+ i^-7r+7r+), and B{D*+ D+j/n^), are prop- 
agated in the fit by adding Gaussian constraints to the 
corresponding fit parameters. The resulting uncertain- 
ties of the measured branching fraction ratios are ex- 
tracted by subtracting in quadrature the statistical un- 
certainties of the fits with branching fraction constrained 
and the one where they are fixed to the central val- 
ues. When calculating the absolute branching fractions 
6(5° ni*^'^ Di*^~) an additional relative uncertainty 
of 16 % is introduced by the measurement uncertainties 
of fs/ fd and the branching fraction of the normalization 
channel 5° — )■ DfD^. The systematic uncertainties are 
summarized in Table [H 

As a result we obtain Jd^d, = 0.183 ± 0.021 ± 0.017, 
f^.^j^ = 0.424 ± 0.046 ± 0.035, /d._d; = 0.654 ± 0.072 ± 
0.065, and fjjMjjM = 1-261 ± 0.095 ± 0.112, where 
the first uncertainties are statistical and the second sys- 
tematic. Taking B(B° D+D-) = (7.2 ± 0.8) x 
10-3 from Ref. [H and fs/fd = 0.269 ± 0.033 from 
Ref. [l3, HO] an absolute inclusive branching ratio of 
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B{B° Di*'>+Di*^~) = (3.38 ± 0.25 ± 0.30 ± 0.56) % is 
calculated where the third uncertainty comes from the 
normalization. Assuming Eq. ([T|) to hold this would 
translate into a decay width difference contribution of 
the B° Di*^+Di*^~ modes of Ar,/r, = (6.99 ± 0.54 ± 
0.64 ± 1.20)%, which is consistent with the standard 
model expectation [2lj . 

In summary, we have measured the branching ratios 
of BO ^ DfDj, BJ^ ^ Dt^Df, 5° ^ D*+D*-, 
and i?" Di*^^Di*^ decays relative to the normal- 
ization mode i?*^ — ^ DfD~. Compared to previous anal- 
yses, we have reduced the systematic uncertainties by 
taking into account the full Df K^K^n^ Dalitz 
structure, as opposed to using a simple two-body D'^ 
decay model. The derived absolute branching ratios 
of 6(Bj ^ DtDj) = (0.49 ± 0.06 ± 0.05 ± 0.08)%, 
6(5^ ^ Dl^DJ) = (1.13 ± 0.12 ± 0.09 ± 0.19)%, 
B{B° D*+D*-) = (1.75 ± 0.19 ± 0.17 ± 0.29) %, and 
i3(B0 -> Di*'^^Di*'^') = (3.38 ±0.25 ±0.30 ±0.56)%, 
where the uncertainties are statistical, systematic, and 
due to the normalization, are the most precise measure- 
ments to date. The central values are lower than but con- 
sistent with the Belle result [ll| and the previous CDF 
result, which is superseded by this measurement. 
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